Excess energy of an ultracold Fermi gas in a trapped geometry 
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We have analytically explored finite size and interparticle interaction corrections to the average 
energy of a harmonically trapped Fermi gas below and above the Fermi temperature, and have 
obtained a better fitting for the excess energy reported by DeMarco and Jin [Science 285, 1703 
(1999)]. We have presented a perturbative calculation within a mean field approximation. 
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I. INTRODUCTION 

Observation of quantum degeneracy of ^"K atoms [l| in 
harmonically trapped geometry made the study of ultra- 
cold fermions an huntingground to the experimentalists 
[ll-Q and theoreticians Q within the last ten years. In 
the remarkable experiment, DeMarco and Jin measured 
the energy of harmonically trapped weakly interacting 
""^K atoms for different temperatures (T) They plot- 
ted the deviation of the average energy of a single 
atom from its classical part (SfcT). This deviation (ex- 
cess energy) particularly below the Fermi temperature 
[Tp ~ 0.6 /^K) clearly indicates that "^^K atoms obey 
Fermi-Dirac statistics. They also compared their experi- 
mental data with the ideal gas prediction. Although this 
prediction matches well with the experimental data yet 
the measured energy is a little higher, and it maintains a 
systematic narrow gap with it in particular for T /Tp < 1. 

It is obvious that confinement in a smaller region and 
repulsive interparticle interaction increase the average 
energy of a particle. Since the thermodynamic limit is 
not well satisfied for our system, as its size (~ mm) and 
number (~ 10^) of particles are small [H, finite size ef- 
fect may contribute in removing the narrow gap between 
the theory and experiment. DeMarco and Jin produced a 
nearly ideal Fermi gas [ll]. But, interparticle interaction 
can not absolutely be neglected in their experiment, and 
it may also contribute in removing the gap 0. 

Besides the study of strongly correlated electron 
(Fermi) gas, weak interaction effects on Fermi gas have 
also been studied for a long time @. Finite size effect 
on an ideal Fermi gas was also studied a few decades ago 
[lOj . While the Fermi gas for those cases are essentially 
homogeneous, the Fermi gas of our interest is inhomo- 
geneous being trapped by an inhomogeneous magnetic 
field. Finite size effect was explored for such a system by 
the authors of Ref. |llll in a way Grossmann and Holthaus 
did for Bose gas They considered the finite size ef- 
fect on our system by adding the zero point energy of the 
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oscillations and a surface contribution (resulted from the 
error of considering discrete energy levels continuum) to 
the bulk density of states [llj. They eventually showed 
that these two corrections cancel each other in the energy 
expression to the first order in angular frequency. Conse- 
quently, they did not get considerable finite size correc- 
tion at all. However, apart from these, finite size effect 
may also be there in measuring temperature p^ . and in 
calculating the chemical potential as well. Incorporating 
these features one can get considerable finite size correc- 
tion to the first order in angular frequency. Along with 
the finite size effect we also consider the contribution of 
weak interparticle interaction in energy for entire range 
of temperature. Although the interest of exploring the 
thermodynamic behavior of an ultracold Fermi gas has 
been moved to strongly correlated regime p^ - [l9j yet the 
same for weakly interacting regime has not surprisingly 
been studied except that for T — > [13, HJli and for 
T 7^ with a particular interest in phase separation of 
its multi-components j23 [. 

In this paper we will incorporate the temperature de- 
pendence of finite size and weak interaction effects on a 
harmonically trapped ultracold Fermi gas in a perturba- 
tive manner. Our calculation will begin in section-II with 
a mean field Hamiltonian where we will consider short 
ranged interaction among the particles [1^ [2l| . We will 
write the average occupation number for the interacting 
fermions in a self consisted way Giorgini, Pitaevskii and 
Stringari wrote for an interacting Bose gas [2^ . We will 
obtain a considerable shift in chemical potential due to 
the interparticle interaction. In section-Ill we will in- 
corporate finite size effect by considering shifts in tem- 
perature and chemical potential apart from the previous 
treatment 11] . In section-IV we will evaluate the interac- 
tion energy. We will get the excess energy in a closed form 
of fugacity (or chemical potential) in section-V. Then 
we will obtain an approximate chemical potential in a 
closed form of temperature. This approximate chemical 
potential will make the excess energy in a closed form of 
temperature. Then we will replot the experimental data 
for evaporation ramp indicating the temperature depen- 
dence of the number of particles [l[ , and will fit the data 
points. We will eventually plot this excess energy for spin 
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t and l particles by considering the fitting formula for 
the evaporation ramp, and will compare our analytic re- 
sult with the experimental data [l[ . In section- VI we will 
compare the interaction energy with the oscillator energy. 
As a corollary, we will obtain a considerable shift in Fermi 
energy due to the interparticle interaction. Finally, we 
will compare the interaction and finite size corrections. 



where z„ = e^^'/^"^ is the fugacity and /j(x) = x — |i -|- 

|j — ... is a Fermi function (or integral). In mathematics 
literature it is known as polylog function (— L?j(-x)) of 
order j. The chemical potential (or fugacity) is to be 
obtained from the constraint that 



(r)d^r. 



(4) 



II. MEAN FIELD THEORY 

Let us consider a weakly interacting trapped Fermi gas 
of N number of 3-D anisotropic harmonic oscillators each 
having spin 1/2. Each single particle state {p, r} is de- 
scribed by a specific momentum (p) and a specific po- 
sition (r) from the center of the trap, and each state of 
course has 2 spin degeneracy (ct =t, \.) . We consider un- 
equal number of particles in the two spin states. It is con- 
venient to begin with an isotropic trap. We will consider 
the anisotropic trap later. Expectation of the effective 
Hamiltonian operator for our system (in the mean field 
level) is [20] 



In an ideal situation {g — 0), no.(p,r) and z^- in Eqn.(2) 
become na\Tp, r) 
Eqn.(4) becomes 



become n'a\Tp,r) and z'^'' respectively. Thus for g — ^ 



(5) 



Shift in chemical potential due to the 
interaction 
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ep°r"-(T(p,r) 



d^pd^r 
(27r;i)3 



Let the shift in chemical potential due to interaction be 
given by (5/1^'' = ^J'a — ■ This shift can now be obtained 
by the Taylor expansion of na{r) in Eqn.(3) about 5 = 



ncr(r)ncr' (r)d'^r,(l) by considering 5za\= e 



) an implicit function of 



where ep,r = ^ + ^mcj^r^ 



m is the mass of each parti- 
cle, w is the angular frequency of oscillations, g = "^^^"^ 
is the coupling constant for the interparticle interac- 
tion {g5^{Y — r') [13, nil), a-s is the s-wave scattering 
length, cr' is the complementary of spin cr, no-(r) = 
/ ^^(p, r)d'^p/(27r?i)^ is the averaged number density of 
spin a particles in equilibrium, and ncr(p,r) is the av- 
erage number of spin cr particles in equilibrium at the 
single particle state {p,r}. This number according to 
the Fermi-Dirac statistics is given by 
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(2) 



where /ig- is the chemical potential for the weakly inter- 



acting fermions, ep., 



2m 



-I- Veff{v) is the mean field 



energy per particle, Ve//(r) — ^muj^r^ + gn^>{r) is the 
effective mean field potential for a spin cr particle [2^ . It 
is to be mentioned that zero point energy can not alter 
the number distribution of particles as any constant shift 
in Ep^r can always be absorbed by fia- It is also to be 
mentioned that we can write Eqn.(2) only in the ther- 
modynamic limit: N^r oo, w — & N„uj^ = constant, 
where stands for total averaged number of particle 
having spin cr. In this limit, any integral over the dis- 
crete levels {n} (of the oscillator energy: {n + 3/2)huj) is 
the same as that over the phase space ({p, r}). Integrat- 
ing n„{p, r) in Eqn.(2) we get 
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na{r) = ^/3/2(z^e 



^-V,Sf{v)/kT 



(3) 



g. Now, to the first order in g and S^Lo'^ we can write 

-(0)/ N 



n^(r) 



where ni°''(r) is the number density for g = Q. For a 

given N^^ integrations of na{v) and /^^'''(r) in Eqn.(6) 
over r are the same. Thus we get 



5(0) 



(r) + 



Sz, 



(i) 



/J(J /\rp 



(6) 



gny{r) 



")d^r = 0. (7) 



Now, to the first order in g and 5za^ this shift is given 
by 



5z^ 
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(-1) 



i+j 



kTXy2iz,) (^+j)3/2i3/2^-l/2- 



(8) 



From Eqn.(8) we can also say that chemical potential 
(or Fermi energy as well) decreases due to the repulsive 
interparticle interaction among the particles. This is not 
a surprize as because increase of energy due to repulsive 
interaction does not necessarily mean chemical potential 
to be increased. We will give qualitative argument in this 
regard in section- VI. However, it should be mentioned 
that the Taylor expansion in Eqn.(6) is applicable only 

1 /3 

for weakly interacting regime {asfij (0) ^1) and not 

1 /3 

for strongly interacting regime {asfiJ (0) > 1). 
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In this section we have highhghted the effect of weak 
interaction not only within mean field approximation but 
also within local density approximation which is embed- 
ded in effective mean field potential [1^ . We will evaluate 
the interaction energy per particle later in section- IV. All 
the calculations in this section has also been done within 
semiclassical approximation due to the fact that we have 
integrated over phase space (from the beginning) instead 
of summing over discrete single particle levels [23] ■ This 
semiclassical approximation is valid only in the thermo- 
dynamic limit. As our system does not approach the 
thermodynamic limit, there must be an error of replace- 
ment of the sum by the integration. In the following we 
will incorporate this error as a finite size correction. We 
will also discuss other features of finite size correction as 
well. 



III. FINITE SIZE EFFECT 



as well. Chemical potential on the other hand, is not 
usually measured (28|, but obtained from Eqn.(5) which 
is also true in the thermodynamic limit. Now we pro- 
pose that actual temperature T^^ and chemical poten- 
tial (/Xct^ ) of a finite system (like ours) are different from 
the measured and calculated value respectively, and they 
must be the same in the thermodynamic limit. 

Interparticle interaction may also change the momen- 
tum distribution, and we are considering its effect as 
a shift in fugacity and not as a shift in temperature. 
Thus apart from a finite size correction part {Sfia ), 
fia"^ also has an interaction correction part {6fj,a^) which 
has already been evaluated in Eqn.(8). We expect that 

(f) — y(f) _y would be positive for smaller confinement 

(i.e. for tighter trap), and 5fia^ = — /^cr would be neg- 
ative for a given No-. Now for g — 0, Eqn.(4) along with 
the integration of ncr(p,r) (in Eqn.(2)) over the surface 

part of the density of states (| (2x^)^ ) becomes [H, 



Apart from the replacement of the summation over the 
discrete levels ({n}) by the integration over ({n}) or that 
over the phase space ({p,r}), finite size correction (to 
total energy) may also arise from the zero point energy of 
all the oscillators {J2rr=f i N^^hw) which should now be 
included in Eqn.(l) [ll|. This replacement of summation, 
however, would impose a correction in density of states, 
such that a surface part (by name) | up to the first 

order in w is to be added with the bulk part (^^fip")- "^^^ 
bulk and surface parts of course correspond respectively 
to the first and second terms of the degeneracy (n^/2 -|- 
in/2+ 1) of the nth level of single particle energy e„ = 
(n -I- |)Sw [1^. These are the two common origins of 
obtaining finite size correction. In the following we are 
going to introduce two more new origins. 



A. Corrections in temperature and chemical 
potential 

Macrostate of our system in thermodynamic equilib- 
rium is usually described by the set {w, T, /io.}, where T 
and /icr (according to the zeroth law) are characteristics 
of an heat and particle reservoir. For our finite system, 
temperature is not measured directly 13]. It is measured 
from the momentum distribution of atoms obtained af- 
ter free expansion switching the magnetic trap off [isj . 
The measured distribution is fitted with the momentum 
distribution formula n<,(p) = (^aU)-^ hnizaG'^" ''^"'^'^) 
which can be obtained by integrating nCT(p, r) over r p7l |. 
This fitting formula is true only in the thermodynamic 
limit, and in the classical limit, its width is proportional 
to -s/T. Effect of zero point motion, however, has not 
been included in this fitting formula. Inclusion of this 
would certainly increase the width and the temperature 



2\ huj 



where Za'^'^ 



AO) 



Sza is the actual fugacity for g = 0. 



Let us further assume that T^^ = T only at the absolute 
zero as for T*^^) — ^ 0, experimentalists may note that the 
width of momentum distribution is the minimum, and 
may set the measured temperature as T = 0. Comparing 
Eqns.(5) and (9) we get the finite size correction in chem- 
ical potential for T = 0, as Sf/J'' = —^huj. Henceforth 

we assume that S^a is the same in all temperature as 
commonly assumed for Bose gas particularly below the 
condensation point (Tc) 12:]. This of course implies that 

Sz^f> = zi^'°^ - 4°^ = -z^°hhuj/2kT. Thus keeping N„ 
in Eqn. (9) fixed, we get the shift in temperature in terms 

of Sza'' as 



T 



' 3z^J'^ fsiza) 



, fto; \ 1 ftw folZr,) , 



For Bose gas, below Tc, 6zjp as well as ST^^^ has been 
assumed to be zero for obtaining a sharp condensation 
point and a considerable shift in Tc [13, [IB]- But, in re- 
ality, a finite system of Bose gas does not have a sharp 
condensation point 12911. Considering the same for Fermi 
gas, authors of Ref.[ll[ also did not get any consider- 
able finite size correction in energy. That is why we have 
assumed Sz^ as well as 6T^^'^ to be nonzero. Our as- 
sumptions do not give a sharp condensation point of a 
Bose gas, and on the other hand, it may considerably 
improve the ideal gas prediction for our ultra-cold Fermi 
gas [H. 
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B. Correction in energy with anisotropy 

Let us now obtain the total energy of the system. In 
an ideal situation {g — 0), Eqn.(l) along with the zero 
point energy and the integration over the surface part of 
the density of states gives the energy of spin a particles 
up to the first order in uj as 



Thus for weak interaction and small number fluctuation, 
we give an ansatz: 



Am) 



m 



So far we have considered our system to be isotropic 
{ujx = ojy = uJz = uj), and we already have noticed 
that tighter trap (larger w) would give a larger finite size 
correction. Finite size corrections for different values of 
u]x,0Jy, UJz would be larger than that for uj^ — — i^z 
due to the fact that anisotropy results smaller (tighter) 
confinement and may even reduce the dimensionality eg. 
for Wz — oo our 3-D trap may become a 2-D one. How- 
ever, the prescription of including anisotropy (to the first 
order) is to replace uj in the finite size correction parts 
by its arithmetic mean w, and in the bulk parts by its 
geometric mean (a)) (soj . Thus we recast Eqn.(ll) along 
with Eqns.(lO) and (5), and get energy per particle in 
terms of measured temperature as 



/3(4°)) 



to / /3(4°V3(4^0 



fc^V-/4(.W)/2(4°^) 2 



(12) 



The inclusion of finite size effects along with the in- 
teraction effects have be done as separate corrections as 
all of their effects in our case are perturbative with re- 
spect to the unperturbed part (i.e. first term) of the 
expectation of Hamiltonian operator in Eqn.(l). Finite 
size effect can not, of course, be treated separately for 
strongly interacting regime. 



-jf^ SO that the interaction 



energy can be evaluated in a closed form. This ansatz, of 
course, is not applicable for strongly interacting regime 
and large number fiuctuation. Now, with this ansatz, we 
get interaction energy per spin a particle as 



(11) 



where 



h{z(J>)) 



V/3(.i°^) 



[hi 



)]Vdx. 



(13) 



(14) 



One check that for T > Tp, the integral representation 
of h{z^a^) in Eqn.(14) can be approximated by the closed 
form 2^[/9/4(4°^)]^ 



V. CORRECTED EXCESS ENERGY 

Energy per particle {Ua) to the first order in ui and 
g can be obtained as a summation of its two parts in 
Eqn.(12) and (13). It is easy to check that its classical 
part is ^ 3kT. Thus we get excess energy per particle 

for our anisotropic system as 

hU"') /3(4°^) 









h{z^') \ 




/3(4°)) y 








-2} + 3 



where kTp gets the value hhj{&N„Y/^ | 



ujfjrpr^\( 

Cj\{m^Y/^)\ /4(z(°)) /,(z(°)) 

- 1, (15) 



IV. EVALUATION OF INTERACTION ENERGY 

Interaction energy per particle up to the first order in 
g can be obtained from Eqn.(l) and (5). It is of course 
difficult to be evaluated for the presence of different z^^ 
and z'f) in fSa\Y) and n^^){v) respectively. At this stage 

fi(0)(r) 

we need to evaluate -fk — specially for weak interaction 

and smaU number fluctuation ( ^"j^^"' = < 1)- 

While for T/Tp — ?> oo, we obviously have the classi- 

cal result: jL -> for T/Tp 0, we have 

^ - 0{^^)) within Thomas- 

Fermi approximation ^8i]. Although for small number 
fluctuation, these two limits do not differ (~ ap- 
preciably, yet we stick with the Thomas-Fermi result as 
because weak interaction plays role only for T/Tp ^ 1. 



A. Chemical potential in a closed form of 
temperature 

If we want to plot the excess energy (in Eqn.(15)) with 
respect to T, we need to know the temperature depen- 
dence of ^^a^ from Eqn.(5). But, it is impossible to be 
known as the inverse of a polylog function does not ex- 
ist. There is no other way but to get /ii"'' in a closed 
form without an approximation. However, whether an 
approximate /ii"-* is good can be justified from the exact 
graphical solutions which might have been used by De- 
Marco and Jin for plotting the ideal gas prediction [l|. 
In the following we prescribe an approximate method for 
obtaining a good ^^a^ in a closed form relevant for the 
entire range of temperature. 

Let us first investigate the approximate /ii^'' {T) in the 
degenerate regime [T/Tp <^ 1). For ^ — ^ 0, Fermi 
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FIG. 1: Dashed and dotted lines represent Sommerfeld's re- 
sult (Eqn. (17)) and improved classical result (Eqn.(19)) re- 
spectively. Solid line represents the approximate chemical po- 
tential in Eqn. (20). Points represent exact graphical solutions 
of yUa(T) in Eqn. (5). 



functions follow Sommerfeld's asymptotic formula [3l| 



r(j + 1) 



1 



6 [Aii°VfcT]2 



(16) 



Now for j = 3, Eqns.(5) and (16) lead to the asymptotic 
form of chemical potential 



1 



TT 



T 



(17) 



One can check that the approximate chemical potential 
in Eqn. (17) is reasonably good for T jTp < 0.5. 

Let us also investigate the approximate /ii"^ (T) in the 
classical regime {T/Tp ^1). Taking only the first term 

of the Taylor series expansion of fsizj^^) in Eqn. (5) we 
get the approximate chemical potential as 



,(0) 



T 



{Tyy)=~3kTF—ln(6^/' — 



Tf 

,(0) 



T 



(18) 



Although the approximate fj,]j' in Eqn. (18) is valid for 
T/Tp — >• oo yet one can check that it is reasonably good 
even for T/Tp > 1.5. 

However, in DeMarco and Jin's experiment [1], ob- 
servation of quantum degeneracy was highlighted for 
0-5 ^ T/Tp < 1.5. In this regime none of the above ap- 
proximations is good. In the next, we are going present 
a calculation of obtaining a more accurate form of the 
chemical potential relevant not only for 0.5 < T/Tp < 1.5 
but for the entire range of temperature. 

Classical result in Eqn. (18) can be improved by consid- 
ering the first three terms of the Taylor series expansion 



of fz{z^^)- Real root for the cubic Eqn. of z'^' now leads 



(0) 



to 



T 

kTp — In 

Tp 



9 165 X 32/3 



+ 



8/W 



(19) 
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T 



FIG. 2: A few experimental data points having 50% uncer- 
tainty are adapted from Ref. [J] . Solid line represents a possi- 
ble fitting (Eqn. (21)) for the data points. 



where f{t) = i/ [128 - 783^^ + 16 V64 - 783^3 -t- 8244i6] i/3 
and t = T/Tp. Sommerfeld's result (Eqn.(17)) and 
the improved classical result (Eqn. (19)) intersects at 
t « 0.468. Thus approximate chemical potential for 
Q < -S- < 00 can be written from the combination of 

— J-F 

these two as 



^i'^\T) = ^i^^'{T<)9{0A68 -t)+ ^i^^nT>)0{t ~ 0.468). (20) 

Comparing with the exact graphical solutions in FIG. 1 
we can say that approximate ^Ij ^ in Eqn. (20) is reason- 
ably good for the entire range of temperature. 



B. Data fitting for evaporation of particles 

However, for comparing our analytic result with the 
experimental data, we need to know the temperature de- 
pendence of the number of particles that arises due to 
evaporation of particles during the experiment of De- 
Marco and Jin [1]. Although the average number of 
particles is not a constant yet they compared their ex- 
perimental data with the ideal gas prediction. Ideal gas 
prediction for a fixed N^- remains unaltered if N^- changes 
quasi-statistically with temperature. Thus temperature 
dependence of N„ in our theory would come only in the 
corrections parts in Eqn. (15). 

DeMarco and Jin showed the temperature dependence 
of Nf in a figure keeping : = 60% : 40% unchanged 
within 5% error during the evaporation ramp [l[. We 
replot their data in FIG. 2, and a possible fitting for these 
data points can be obtained with the help of Wolfram 
Mathematica, and it is given by 

N^{t) = 10"(CO + CltTS + C2i5 + C3<3 + + c^t), (21) 

where cq = -2.05293, a = 6.55557, ca = -6.03520, = 
1.69586, C4 = -0.168018, cg = 0.00472971. From FIG. 2 
we can see that varies from ^ 2.7 x 10^ to ^ 4.2 x 10^ 
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T 



FIG. 3: Dotted line represents ideal gas prediction Solid 
and dashed lines represent our result in Eqn.(15) for spin f 
and 4- respectively. Several experimental data points having 
20% uncertainty are adapted from Ref. Here we take m — 
40 X 1.675 X 10"^'^ kg, cD/w = 1.3677, Nf : = 3 : 2 and 
as = 169ao [3. 



as T/Tf decreases from 2.7 to 0.5. Multiplying the right 
hand side of Eqn.(21) by 2/3 we get Ni{t) as well. 

Apart from the explicit temperature dependence of the 
different parts in Eqn.(15), finite size correction part 
in this equation also has an implicit temperature de- 
pendence as DeMarco and Jin smoothly varied the ra- 
dial trap frequency cjj_ {— lu^ — cuy) from 27r x 44 Hz 
to 27r X 370 Hz keeping the axial trap frequency fixed 
{ujz — 2tt X 19.5 Hz) for obtaining the temperature de- 
pendence of Ncr ilj. Thus ujj_ in uj/uj of Eqn.(15) has an 
implicit temperature dependence which was not shown in 
Ref.[l|. However, for a quasi-statically smooth variation 
of one can expect the ratio of arithmetic and geo- 
metric means of ujj. , ujy and ojj, to be a constant around 
T/Tp = 1, and we set uj/Cj = 1.3677 by picking up 
= 27r X 137 Hz for T/Tp « 1 [J. 



C. Our result for excess energy per particle 



We replace m different parts of Eqn.(15) by Na-{t) 
in Eqn.(21), and plot the excess energy in Eqn.(15) 
(along with the approximate bulk chemical potential in 
Eqn.(20)) in FIG. 3 for spin f and spin I particles keep- 
ing uj/uj unchanged with its value at T w Tp. In the 
same figure we also plot the ideal gas prediction putting 
TVct = oo and g = in Eqn.(15). From this figure we 
can say that our analytic result matches better with the 
experimental data. 



VI. COMPARISONS AMONG IDEAL AND 
CORRECTION PARTS 

As a consistency check and for showing plausibility of 
our theory, let us now compare the interaction energy 
and potential energy for different temperatures. It is easy 
to understand from Eqn.(13) that N„Ua'^^'' is the total 
interaction energy We can also say from the 

virial theorem that {N„u'^^ + N„.u'^})/2 is the total 
oscillator energy {E'^^°^). Thus from Eqn.(15), we get 
their ratio for N„ = N^' = N/2 as 



£;(ho) 



9 



(22) 



We plot the right hand side of Eqn.(22) for a fixed num- 
ber of particles in FIG. 4. However, to obtain ^(T^ for 
T 0, we must need to know the asymptotic expan- 
sions of f4{z^^) and h{z^a ^). While the first one is fol- 
lowed from Eqn.(16) second one is to be obtained from 
Eqn.(14). In the integral representation of h{z^'^), the 

Fermi function contributes essentially for x < (Iuzct"'')^^^. 
Now from Eqns.(14) and (16) we get the desired asymp- 
totic expansion as 



1024(k 



28357r3/2 



1 



2171^ 



32 (Ir 



(23) 



This form has been used in FIG. 4 in particular for 
T/Tp < 0.1. Thus for T ^ 0, Eqn.(22) gives 

£;(int)(o) _ 8192 X \/2 X 31/6 /ivi/6a 



£;(ho)(0) 



0.497 



28357r2 



I 



(24) 



where I = y is the confining length scale of a particle 

in our anisotropic trap. One can check that Eqn.(24) is 
consisted with the Thomas- Fermi result for T — ^ :20] . 

As a corollary of the asymptotic expansion, we can also 
calculate the shift in Fermi energy (from Eqn.(8)) due to 
interaction as 



128 X 72 X 3^/^ {N^/^a, 



-0.210^ 



1057r2 



I 



(25) 



That Fermi energy (or Fermi temperature as well) de- 
creases in Eqn.(25) is not a surprise as because the num- 
ber density of fermions at the central region of the trap 
decreases for repulsive interaction. It is a similar phe- 
nomenon of how the condensation temperature of a har- 
monically trapped Bose gas decreases for repulsive inter- 
action [30|. It should also be mentioned in this regard 
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FIG. 4: Dotted line represents the ratio of interaction and 
potential energy (Eqn.(22)) for fixed number of particles in 
units of UsN^^^ /I. Solid and dashed lines represent the ratio 
of interaction and finite size corrections (in Eqn.(15)) for spin 
t and J, particles respectively. 



that Pauli paramagnetism arises for /i4.(0) > MtC^) while 

[/t(o) > UM UM- 

It is not clear from FIG. 3 how big is the finite size 
correction with respect to the interaction correction. Di- 
viding ?7ct"*^ by C/a^"'' in Eqn.(15) we get the ratio of 
the interaction and finite size correction, and plot it in 
FIG. 4 for spin f and \. particles with the experimen- 
tal parameters used in FIG 3. Now we can say, from 
FIG. 4, that interaction and finite size corrections are 
comparable within the temperature range of our interest. 
This comparison for T /Tp < 0.5, however, is ill defined 
as the temperature dependence of is insensitive for 
T/Tf < 0.5 in FIG. 2. 



approximation as a better choice with respect to our 
choice for the mean field (Hartree-Fock) approximation 
• However, we can not get any contribution from the 
ladder diagrams if we consider short ranged interaction 
{gS^{r — r')) among the particles. Hence, Hartree-Fock 
approximation is reasonably good even for the low dense 
gas considered by us. 

Finite size correction scheme prescribed here can also 
be generalized for the Bose gas. For Bose gas, Fermi 
functions (/j(x)) are to be replaced by the Bose-Einstein 
functions {gj{x)). 

Approximate chemical potential (in a closed form) ob- 
tained by us makes an analytic theory plausible in the 
entire range of temperature. 

Our main new result is the finite size correction part 
(which is first order in (|^) ~ 0.01) in Eqn.(15). It has 
come from the considerations of (i) zero point energy, 
(ii) discreteness of single particle energy spectrum, (iii) 
shift in temperature, and (iv) shift in chemical potential. 
While first two considerations result correction to the sec- 
ond order in (^) the last two considerations result 
to the first order in . This first order correction along 
with interparticle interaction in Eqn.(15) may consider- 
ably improve the ideal gas prediction. Our result can, 
of course, be negligibly improved to the second order by 
including the result of Ref. [ll| . Instead of taking an in- 
termediate value of (D/w, our result can also be improved 
a little if one gets the data for the implicit temperature 
dependence of uj±_ (or of uj/lj) in Eqn.(15). 

Excess energy in other experiments also found to be a 
little larger with respect to the ideal gas prediction 
Our result (for the finite size and interaction corrections) 
are somewhat uncertain at the present level of accuracy 
[l], [2] , and compelling the experimental evidences for its 
presence. Future experiment with reduced error may con- 
firm our prediction. 



VII. CONCLUSION 



We have presented an analytic calculation for explor- 
ing thermodynamic behavior of a harmonically trapped 
weakly interacting Fermi gas through the temperature 
dependent energy. Our result for ^ — >■ and ^ 3> 1, 
reproduces the standard Thomas-Fermi and classical 
results [H respectively, and of course matches better with 
the experimental data 

Since the number density of the fermions in the trap 
is reasonably low (~ 10^^/m^) one may expect ladder 
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